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Abstract. We discuss the different definitions of the mass 
of a halo in common use and how one may convert between 
them. Using N-body simulations we show that mass esti- 
mates based on spherical averages are much more tightly 
correlated with each other than with masses based on the 
number of particles in a halo. The mass functions pertain- 
ing to some different mass definitions arc estimat ed an d 
compared to the 'universal form' of Jenkins et al. ( |2000|) . 
Using a different simulation pipeline and a different cos- 
mological model we s how t hat the mass function is well fit 
by the Jenkins et al. (200C) fitting function, strengthening 
the claim to universality made by those authors. We show 
that care must be taken to match the definitions of mass 
when using large N-body simulations to bootstrap scal- 
ing relations from smaller hydrodynamical runs to avoid 
observationally significant bias in the predictions for abun- 
dances of objects. 
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1. Introduction 

One of the most fundamental predictions of a theory of 
structure formation is the number density of objects of 
a given mass, the mass function, at a given redshift. Ac- 
curate mass functions are used in a number of areas in 
cosmology; in studies of galaxy formation, in measures of 
volumes (e.g. galaxy lensing) and in attempts to infer the 
normalization of the power spectrum and the density pa- 
rameter from the abundance of rich clusters. In the latter 
case the mass function is the point of contact allowing 
us to bootstrap the excellent statistics of large N-body 
simulations with observable properties of clusters normal- 
ized for example by hydrodynamic simulations of smaller 
volumes. In this way we can obtain reliable estimates of 
e.g. the number of clusters as a function of temperature 
and redshift, which can in turn be used to constrain the 
matter density, the normalization of the power spectrum 
and the statistics of the initial density field. 

Several different definitions for the "mass" of a halo 
are in common use, each having different advantages. For 



example there are at least 3 different definitions of mass 
for the mass-temperature relation as computed from hy- 
drodynamical simulations of galaxy clusters and all are 
different from the definition of mass commonly used in 
the mass function which is itself different from the mass 
usually employed in analytic studies based on the Press- 
Schechter ( |1974| ; hereafter PS) theory. Observational data 
have improved to the point where it is important to dis- 
tinguish between these different definitions of mass, lest 
we bias our theoretical predictions. We give an example 
of how such a conversion can be made, at least approxi- 
mately, for a certain class of mass estimator. 

2. Halo mass definitions 

In this section we describe some of the definitions of the 
mass of a dark matter halo in common use. This list is ob- 
viously not exhaustive, but it is representative. We begin 
by recalling some background about the spherical top-hat 
collapse model (see e.g. Peebles 1993; Peacock 199£; Lid- 
die & Lyth 200C , and references therein) from which much 
of the language in this field has been borrowed. 

The spherical top-hat ansatz describes the formation of 
a collapsed object by solving for the evolution of a sphere 
of uniform overdensity 8 in a smooth background of den- 
sity p. By Birkhoff 's theorem the overdense region evolves 
as a positively curved Friedman universe whose expansion 
rate is initially matched to that of the background. The 
overdensity at first expands but because it is overdense the 
expansion slows (relative to the background) and eventu- 
ally halts before the region begins to recollapse. Techni- 
cally the collapse proceeds to a singularity but it is as- 
sumed in a "real" object virialization occurs at twiceQ the 
turn-around time, resulting in a sphere of half the turn- 
around radius. In an Einstein-de Sitter model the over- 
density (relative to the critical density) at virialization is 
A c = 187r 2 ~ 178. We shall always use A c to indicate the 
overdensity relative to critical of a virialized halo, which 
will be lower for smaller Om- Note that some authors use a 
different convention in which A c is specified relative to the 

1 There is a small correction to this in the presence of a 
cosmological constant which contributes a Ar 2 potential. 
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Fig. 1. The relations between various definitions of the 
mass of a halo as a function of Dm assuming the halo 
density profile follows the NFW form with concentration 
parameter c = 5. 



background matter density - our A c is Dm times theirs. 
The linear theory extrapolation of this overdensity is nor- 
mally denoted S c and is (3/20)(12tt) 2 / 3 ~ 1.686 in an 
Einstein-de Sitter model. This overdensity is often used 
as a threshold parameter in PS theory and its extensions 
and has a very weak cosmology dependence which is often 
neglected. 

We now turn to some of the definitions of mass in use 
in the literature. Let us suppose that we have identified 
some object, a problem to which we shall return in §|], 
and have chosen a fiducial 'center' about which to take 
spherically averaged profiles. We define Ma as the mass 
contained within a radius r& inside of which the mean 
interior density is A times the critical density 



r dr p(r) 



A 



"Pcrit^A 



(1) 



The 'virial mass' from the spherical top-hat collapse model 
would then be simply Ma c . We shall refer to this mass as 

Mth-vir- 

An alternative interpretation of the spherical collapse 
model is that the virial radius corresponds to that point 
within which the material is virialized and external to 
which the mass is still collapsing onto the object. Some 
simulations suggest that this occurs at A = 200 more or 
less independent of cosmology, and so a common mass es- 
timator is M200 which is approximately the virial mass if 
Qm = 1- Authors also use M500 or M1000 in some appli- 
cations. 

The final mass estimator we shall consider is the one 



used by Jenkins et al. ( 2000 ) in their analysis of simula- 
tions performed by the Virgo consortium. The mass they 
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Fig. 2. Correlations between two estimators of the mass, 
Aff o f0.2 and M2001 m simulations of a ACDM model (see 
text). Only halos with more than 1000 particle are used. 
The upper panel shows that the different estimates of mass 
are correlated, but offset (the dotted line marks y = x). 
The dashed line indicates the prediction for the offset of 
The lower panel shows the mass ratios as a function of 
M200 on a linear scale. 



assign to a halo is simply the sum of the masses of the 
particles identified as members of the halo by their halo 
finder. This mass estimator is the least easily interpreted 
theoretically, though it is well defined algorithmically, but 
has the very convenient feature that the mass function for 
this estimator is independent of cosmology! 

Unfortunately there is no unique algorithmic defini- 
tion of a dark matter halo, even within a 3D simulation 
itself. To define the mass of an object we must first specify 
the object in question! Although other halo finders are in 
common use, we shall deal exclusively with halos found us- 



ing the Friends-of- Friends (Davis et al. 1985) algorithm, 
hereafter called FOF. Our reason for this choice is that 
it is the group finder employed by the Virgo consortium 
in deriving the mass functions from their very large N- 
body simulations (Jenkins et al. |2000 ) and in terms of 
which they find a universal mass function. The FOF al- 
gorithm has one free parameter, b, the linking length in 
units of the mean inter-particle spacing. Commonly used 
values of b are 0.1, 0.15 and 0.2, although other choices 
exist (e.g. Gardner 2000 chooses instead b~ 3 = f2]y[A c /3). 



Jenkins et al. (2000) find that the mass function is uni- 
versal if they take b = 0.2, indepe ndent of the cosmology 
under consideration. Jenkins et al. ( 2000 ) claim that in the 
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MfofO.2 


A/th-vir 


A/200 


MfofO.2 


1.000 


0.974 


0.966 


Afth-vir 


0.974 


1.000 


0.997 


M 200 


0.966 


0.997 


1.000 



Table 1. Correlations between the various mass defini- 
tions for the halos of our large simulation. Only halos with 
more than 1000 particles are used. 



limit of a large number of particles per halo, FOF finds 
all particles within an iso-density contour b~ 3 times the 
mean matter density or ^m/^ 3 times critical density. We 
shall return to this point in §|[ 

In the next sections we describe how these different 
masses are related. 



3. An approximate conversion 

It would be advantageous if we could convert, at least 
approximately, between the different definitions of mass 
in common use. In this section we make a first attempt 
in this direction. We assume that all halos are spherical 
and have a universal profi le, fo r example^ the NFW form 
(Navarro, Frenk & White |l996| ): 



p(r) (xx 1 (1 + x) 



(2) 



where x = r/r s and r s is a scale radius usually specified in 
terms of the concentration parameter c = r^ooAs- Navarro 
et al. ( 1996 ) refer to r2oo and Ahoo throughout as the 
'virial radius' and 'virial mass' respectively. Again, N-body 
simulations have shown that the concentration parameter 
is a weak function of virial mass, taking the value c ~ 5 for 
masses characteristic of clusters which shall be our focus 
here. 

With a given profile it is straightforward to relate the 
various mass definitions as shown in Fig. |l|. Note that 
the masses can differ significantly in this model and that 
the cosmology dependence is quite strong. How well does 
this crude model predict the relationship between different 
mass definitions in practice? 

4. Simulations 

The model of the previous section made several approxi- 
mations. To compare these different mass estimates with 
each other more carefully we have used some N-body sim- 
ulations, originally run for another purpose. We simulated 
the Ostriker & Steinhardt (1995) concordance model, 
which has Q m = 0.3, fl A = 0.7, H = 100 hkms^Mpc^ 1 



2 There is some controversy in the literature about the slope 
of the density profile as r — > 0. Luckily departures from the 
NFW form are expected to occur only at very small radii, con- 
taining a vanishingly small amount of the total mass, and thus 
do not concern us here. 



with h = 0.67, n B = 0.04, n = 1 and cr 8 = 0.9 (cor- 
responding to Sh = 5.02 x 10~ 5 ). For this cosmology 
A c ~ 101 from the top-hat collapse model. 

We have used two 'high' resolution and two lower res- 
olution N-body simulations. The first was a 256 3 particle, 
dark matter only simulation in a 100 3 ft. _1 Mpc box. The 
simulation was run on 16 processors of the Origin2000 at 



NCSA with the TreePM-SPH code (White et al. |2000| ) op- 
erating in collisionless (dark matter only) mode. The grav- 
itational force softening was of a spline form (e.g. Hern- 
quist & Katz 1989| ), with a "Plummer-equivalent" soft- 
ening length of 15/i _1 kpc comoving. The simulation was 
evolved from z = 100 until the present and took a to- 
tal of 1300 (normalized) CPU hours. A second 256 3 par- 
ticle simulation was run from z — 70 in a 256/i~ 1 Mpc 
box with a 35/i _1 kpc softening. This run took a total of 
880 (normalized) CPU hours, also on 16 processors. To 
provide additional statistics on the high-mass end of the 
mass function we additionally ran several smaller (128 3 
particles) simulations of the same model in boxes of side 
200/i~ 1 Mpc. 

From the z = outputs we generated group catalogues 
using the FOF algorithm with b = 0.2. For each halo we 
defined the center for our spherical averages as the parti- 
cle with the minimum potential energy. This corresponds 
closely to the most bound particle and the density peak for 
a halo in all but the most disturbed systems, and is more 
robust than the center of mass. With the group catalog 
and centers so defined it is straightforward to calculate 
each of the estimators described above. Throughout we 
shall use M200 as our "base" mass to which the others are 
compared, since this fits best into the philosophy of §||. 

We show in Fig. || that the different mass estimators 
are indeed highly correlated as one would expect (see also 
Table |lj). The ratio of the FOF mass to M200 is seen to 
span a reasonable range and be offset from unity by a non- 
negligible amount. Fig. |^ shows a histogram of the ratio 
of Mf fo.2 and M t h- V ir to A/200 for all of the halos above 
a given number of particles. 

Let us focus first on M t h- V ir- The shape of the his- 
togram is roughly the same regardless of particle number 
indicating that the distribution is not strongly affected 
by finite particle numbers in the halos chosen. The top- 
hat virial mass and M200 are fairly tightly correlated, as 
one might expect since they are both spherically averaged 
statistics and the density contrasts are not too different. 
The approximation of ^ provides a reasonable estimate of 
the ratio of masses, presumably because of the spherical 
averaging being performed. 

The situation with regards Mf fo.2 is more compli- 
cated. Firstly the scatter is much larger, since the groups 
found by FOF can be quite irregularly shaped (see Fig. [3]) . 
The actual ratio of A/f o f .2 to M 2 oo is larger than predicted 
by our model. By eye the halo finder looks to have merged 
close groups. It is possible that with better mass resolu- 
tion this would occur to a lesser extent, although the lack 
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Fig. 3. Distribution of mass ratios for halos from the first simulation. The solid line shows all halos with more than 
500 particles, the dashed line 1000 particles and the dotted line 2000 particles. The ratios predicted by the model of 
§|are M ioi0 . 2 /M 20Q = 1.07 and M th _ vir /M 200 = 1.22. 



of a clear trend with particle number in Fig. ^ argues that 
convergence would be quite slow. 

To investigate whether FOF has merged adjacent halos 
and whether a FOF halo can be considered as all particles 
above a given density threshold (as assumed in §||) we 
have examined the particles in the two groups of Fig. |]. 
For each particle we calculate its distance from the halo 
center and its density, defined from the distance (r ne i g ) to 
the 32nd nearest neighbour^] as 



3M 

47T 



(3) 



where M is the enclosed mass. This density estimate is 
somewhat noisy, but will be sufficiently accurate for our 
purposes. 

The results are shown in Fig. along with the appro- 
priate NFW profile (assuming c = 5) for comparison. Note 
that the substructure in the halos is readily apparent as 
density 'spikes' as functions of radius and can be matched 
to 'sub-halos' in Fig. || 

It is clear from Fig. || that FOF is (at least approxi- 
mately) linking all particles above a density Om6 -3 of crit- 
ical. However the chosen contrast is low enough that this 
includes particles well outside r 2 oo (or the virial radius) 
and the presence of these 'nearby' particles is biasing the 
masses above M200 and contributing to the extra scatter. 
This problem could be mitigated to some extent by reduc- 
ing b, but it is not clear whether the mass function would 
remain universal if a different (cosmology independent) b 
were chosen. 



Finally we show in Fig. ^ the different mass functions 
that would be derived from these halos using the 3 differ- 
ent mass definitions. Note that our results, obtained with 
a completely independent N-body code and analyzed with 
a completely independent set of software, agree well with 
the universal mass function of Jenkins et al. ( 2000 ). How- 
ever, especially at the high mass end, the number density 
is quite sensitive to the mass definition used. Intriguingly 
the fitting formula quoted by Jenkins et al. ( 2000 ) pro- 
vides a good fit to the mass function if we use M t h- V ir as 
our definition of mass, an even better fit than if we use 

Mf o f0.2- 



5. An example: the temperature function 

An example of where this mass uncertainty may be im- 
portant is in using the observed temperature function of 
local X-ray clusters to constrain the normalization of the 
present day power spectrum. Although Pen ( 1998| ) has 
argued for a direct prediction of the temperature func- 
tion from simulations, most people use a hybrid approach 
in obtaining this constraint. First a mass function, usu- 
ally calib rated by simulations such as those of Jenkins et 
al. ( [2000 ), is assumed and then a mass temperature rela- 
tion from hydrodynamic simulations is used to predict the 
temperature function. 



3 We use all of the particles in the simulation, not just group 
members, in computing r ne ig- 



Hydrodynamic simulations show good agreement for 
the total mass and X-r ay tem perature properties of clus- 
ters (e.g. Frenk et al. 1999 ). The mass-temperature re- 
lation from these simulations generally follows the virial 
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Fig. 5. The halo density profile, as a function of radial 
distance from the group center, for the groups of Fig. ^j. 
Each dot represent a particle in the group (only a small 
fraction of the particles are shown for clarity) and the 
density is estimated through Eq. |[ The dashed horizontal 
line is JIm^ - , the arrows mark r 2 oo and ^th-vir and the 
solid line is the 'model' NFW profile of §0, with c = 5. 



e.g. Henry (|2000[) for a recent compilation) but in addition 
to this different authors choose different conventions for 
M. For reasonable CIm, the difference in mass definitions 
is about the same size as the difference in proportionality 



10- 



Fig. 4. The mass distribution for two groups from the 
first simulation, along with the top-hat 'virial' radius 
(solid) and r2oo (dashed line). An x — y projection, with 
axes in /i _1 Mpc, is shown centered on the point of min- 
imum potential. The upper panel shows a group with 
M 200 ^ lO 14 h _1 M (62,000 particles; 5% shown), while 
the lower panel shows the largest group in the simulation 
(M 200 = 8 x 1O u /i -1 M 0) 2 x 10 5 particles; 2% shown). We 
have plotted only a fraction of the points for clarity, this 
suppresses the appearance of substructure in the group. 
There is a tendency for the larger groups to be more ir- 
regular. 

relation quite well. Ignoring a small correction from the 
Ar 2 potential, for an object virialized at a redshift z 



M(T,z) oc [T/(l + z)] 3/2 



(4) 



and the proportionality constant is then fixed by the sim- 
ulations. There still remains some disagreement over the 
precise coefficient and its dependence on cosmology (see 
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Fig. 6. The cumulative mass function for the Concordance 
model. The points are mass functions from the simulations 
where the mass is defined as: FOF0.2 (squares), Mth- V ir 
(circles) and M200 (triangles). The solid line shows the 
predictions of the Press-Schechter theory, the dashed line 
is the fit of Jenkins et al. (|2000[). 
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constants between groups and the la scatter in M at fixed 
T found by the different groups. Unfortunately this alone 
does not account for all of the discrepancy. 

To take only two relations as examples, Evrard, Met- 
zler & Navarro ( |1996fl define the mass in Eq. ^ as M2oo- 
The M of Bryan & Norman ( p.998[ ), on the other hand, 
is the 'virial' mass in the sense of the spherical top-hat 
model. In both cases there is ample justification for the 
definition chosen, and the authors are explicit about their 
choice. However, there is an J7m dependent ratio between 
these two definitions which cannot be neglected if the tem- 
perature function is to be predicted to the factor of 2 level! 



6. CoAclusions 



ter for Supercomputing Applications, University of Illi- 
nois, Urbana-Champaign. 
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estimated by assuming a 'universal' spherical profile (§g). 
The hydrodynamical simulations which calibrate observ- 
ables as a function of cluster mass typically use such a 
spherically averaged mass definition. 

Unfortunately these definitions are not very tightly 
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